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Abstract. Recently, antihydrogen atoms were trapped at CERN in a magnetic 
minimum (minimum-B) trap formed by superconducting octupole and mirror 
magnet coils. The trapped antiatoms were detected by rapidly turning off these 
magnets, thereby eliminating the magnetic minimum and releasing any antiatoms 
contained in the trap. Once released, these antiatoms quickly hit the trap 
wall, whereupon the positrons and antiprotons in the antiatoms annihilated. 
The antiproton annihilations produce easily detected signals; we used these 
signals to prove that we trapped antihydrogen. However, our technique could 
be confounded by mirror-trapped antiprotons, which would produce seemingly- 
identical annihilation signals upon hitting the trap wall. In this paper, we 
discuss possible sources of mirror-trapped antiprotons and show that antihydrogen 
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and antiprotons can be readily distinguished, often with the aid of applied 
electric fields, by analyzing the annihilation locations and times. We further 
discuss the general properties of antiproton and antihydrogen trajectories in this 
magnetic geometry, and reconstruct the antihydrogen energy distribution from 
the measured annihilation time history. 



PACS numbers: 52.27.Jt, 36.10.-k 52.20.Dq 
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Figure 1. A schematic, cut-away diagram of the antihydrogen production and 
trapping region of the ALPHA apparatus, showing the relative positions of the 
cryogenically cooled Penning-Malmberg trap electrodes, the minimum-B trap 
magnets and the annihilation detector. The trap wall is on the inner radius 
of the electrodes. Not shown is the solenoid, which makes a uniform field in z. 
The components are not drawn to scale. 



1. Introduction 

Recently, antihydrogen (H) atoms were trapped in the ALPHA apparatus at CERN 
[U [2]. The abihty to discriminate between trapped antihydrogen and incidentally 
trapped antiprotons was crucial to proving that antihydrogen was actually trapped 
[U m [3] . The antihydrogen was trapped in a magnetic minimum 4J created by an 
octupole magnet which produced fields of 1.53 T at the trap wall at i?w = 22.28 mm, 
and two mirror coils which produced fields of 1 T at their centers at z = ±138 mm. 
The relative orientation of these coils and the trap boundaries are shown in Figure [l] 
These fields were superimposed on a uniform axial field of 1 T O [6] . The fields thus 
increased from about 1.06 T at the trap center (r = z = mm), to 2 T at the trap 
axial ends (r = 0mm, z = ±138 mm), and to ^1-062 + 1.53^ T — 1.86T on the trap 
wall at (r — Rw, z = mm). |f] Antihydrogen was trapped in this minimum because 
of the interaction of its magnetic moment with the inhomogeneous field. Ground state 
antihydrogen with a properly aligned spin is a low field seeker; as its motion is slow 
enough that its spin does not flip, the antihydrogen is pushed back towards the trap 
center [|] by a force 



where B is the total magnetic field, and /xjj is the antihydrogen magnetic moment. 
Unfortunately, the magnetic moment for ground state antihydrogen is small; the trap 
depth in the ALPHA apparatus is only ^Tiap = 0.54 K, where K is used as an energy 
unit. 

Trapped antihydrogen was identified by quickly turning off the superconducting 
octupole and mirror magnetic field coils. Any antihydrogen present in the trap was 
then released onto the trap walls, where it annihilated. The temporal and spatial 
coordinates of such annihilations were recorded by a vertex imaging particle detector 

f Note that 0.06 T is field from the mirrors at ^ = mm. 

§ Because of the interaction between the mirror and octupole fields, the magnetic field minimum is 
actually slightly radially displaced from the trap center, not at the trap center itself. 



F = V(/Xh-B), 



(1) 
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OlIllH]. The detector is sensitive only to the charged particles produced by antiproton 
annihilations; it cannot detect the gamma rays from positron annihilations. Thus, 
it cannot directly discriminate between antihydrogen and any bare antiprotons that 
might also be trapped. We must use additional means to prove that a candidate 
observation (event) results from an antihydrogen annihilation. 

Bare antiprotons can be trapped by the octupole and mirror fields because they 
may be reflected, or mirrored [9] , by the increasing field as they propagate away from 
the trap center. Antiprotons obey the Lorentz force, 

F = -g(E + vxB), (2) 

where q is the unit charge, v is the antiproton velocity, and E is the electric field, if 
any, present in the trap. In our circumstances, the antiprotons generally satisfy the 
guiding center approximation requirements 'lOi . Temporarily ignoring E, the force 
law for the antiprotons reduces to one similar to that for antihydrogen, ([T]), 

F = V(Mp.B), (3) 

with the antiproton perpendicular magnetic moment fip replacing /xg in (ij), and 
with the additional constraint that the antiprotons follow the magnetic field lines, 
slowly progressing between lines as dictated by the other guiding center drifts. Here, 
Hp = \fjip\ = £±/B, /Xp is aligned antiparallel to B, and E± is the antiproton kinetic 
energy perpendicular to B. Because /ip is adiabatically conserved, antiprotons can be 
trapped if their parallel energy is exhausted as they propagate outwards from the trap 
center. The trapping condition comes from the well-known magnetic mirror equation, 

Smax =bJi+ ^\ , (4) 



^_L0 

which defines the largest total magnetic field Bmax to which an antiproton that starts 
at the trap center can propagate. Here, Bq is the total magnetic field magnitude at 
the trap center, and £±q and £j[o are the antiproton's kinetic energies perpendicular 
and parallel to the total magnetic field at the trap center. Using Q, we can readily 
define the critical antiproton trapping energy ratio 

i>wall — i'O 

where i?waii is the smallest total magnetic field magnitude on the region of the trap 
wall accessible from the trap center[||] An antiproton will be trapped if its ratio of 
fj_o/f||o exceeds TZpc, i-e., if its perpendicular energy is large compared to its parallel 
energy. 

A typical antihydrogen synthesis cycle [1] starts with 15,000-30,000 antiproton^ 
trapped in an electrostatic well, and several million positrons trapped in a nearby 
electrostatic well of opposite curvature. This configuration is called a double well 
Penning-Malmberg trap; the electrostatic wells provide axial confinement, and the 
aforementioned axial magnetic field provides radial confinement. The antiprotons 
come from CERN's Antiproton Decelerator (AD) [12 , and the positrons from a Surko- 
style [13] positron accumulator. The reader is referred elsewhere jUt i3j for the details 
of the trap operation; we will here discuss only those aspects of the trap operation 
relevant to discriminating between antiprotons and antihydrogen. 

II More completely, -Bwall is the lesser of the total magnetic field at the trap wall or the total magnetic 
field on the trap axis directly underneath the mirror. In our case, the former is lower. 
^ The lower number (15,000) characterizes the number of antiprotons when we employ antiproton 
evaporative cooling 
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About one third of the antiprotons convert to antihydrogen on mixing with 
positrons 2J. Some of these antihydrogen atoms hit the trap wah and annihilate. 
Others are ionized by cohisions with the remaining positrons or antiprotons [3], or by 
the strong electric fields present in mixing region [JM I15j , and turn back into bare 
antiprotons (and positrons). Only a very few antiatoms are trapped at the end of 
the mixing cycle, and confined with these few are approximately 10,000-20,000 bare 
antiprotons. If these antiprotons were isotropically distributed in velocity, it is easy 
to show by integrating over the distribution that the fraction that would be trapped 
by the octupole and mirror fields alone once the electrostatic fields are turned off is: 



Since TZpc = 1-35 for our magnet system, 65%, of an isotropically distributed 
population of antiprotons would be trappedj^ The actual distribution of the bare 
antiprotons is unknown and likely not isotropic. Nevertheless, if any fraction of 
these antiprotons were actually present in the trap when the magnets are shut off, 
the signal from these antiprotons would overwhelm the signal from any trapped 
antihydrogen. Thus, our goal is two-fold: (1) in the experiment, to eliminate the 
trapped antiprotons if possible, and (2) in the analysis, to be able to discriminate 
between trapped antihydrogen and any mirror-trapped antiprotons that might have 
survived the elimination procedures. 

In Section [2] of this paper, we describe the numeric simulations that we used to 
investigate these issues. In Section [3j we describe how we apply large electric fields 
which clear all antiprotons with kinetic energy less than about 50 eV. In Section |4] 
we consider the various mechanisms that could result in mirror-trapped antiprotons 
with this much energy and conclude that few, if any, antiprotons are trapped. In 
Section [5] we describe experiments which benchmark the antiproton simulations, and 
in Section [6] we discuss the postulated antiproton energy distribution. Finally, in 
Section [7] we employ simulations to show that if any mirror-trapped antiprotons were 
to survive the clearing processes, they would annihilate with very different temporal 
and spatial characteristics than do minimum-B trapped antihydrogen atoms. 

2. Antiproton and antihydrogen simulations 

In this section, we first describe how we calculate the electric and magnetic fields 
present in the apparatus, including the effects of eddy currents while the magnets are 
being turned off (shutdown). Then we describe the simulation codes that use these 
fields to determine the antiproton and antihydrogen trajectories. 

2.1. Fields 

Electric fields are generated in the trap by imposing different potentials on the 
trap electrodes (see Figure [T]). In the simulations, these fields are determined by 
finite difference methods. Two independent calculations were undertaken. The first, 
and the one used in the majority of the simulations, was hand coded and used 
a slightly simplified model of the electrode mechanical structure; the second was 

+ This calculation assumes that the antiprotons originate at the magnetic minimum in the trap. 
The parameter, T^pc, is greater for antiprotons that originate elsewhere, so, for such antiprotons, the 
trapping fraction would be less. 
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obtained using the COMSOL Multiphysics package [16. and an exact model of the 
electrode mechanical structure. When the calculations were compared, the largest 
differences in the potentials were near the gaps between the electrodes at the trap 
wall. These differences reflected the handling of the computational grid near the 
electrode gaps. The largest potential energy differences were more than two orders of 
magnitude smaller than the antiproton energy scale. Away from the electrode gaps, 
these differences were more than four orders of magnitude smaller. The annihilation 
location statistics that result from the two finite difference calculations agree within 
•\/]V fluctuations. 

Four magnetic field coils, a solenoid, two mirrors and an octupole, produce the 
fields modeled in the simulations. (A fifth coil present in the experiment, a solenoid 
which boosts the magnetic field during the antiproton catching phase, is not energized 
during the times studied in the simulations.) No simple analytic expressions for 
the field from these coils exist because their windings possess an appreciable cross- 
sectional area and are of finite length. Consequently, we use the Biot-Savart numeric 
integrator found in the T0SCA/0PERA3D field solver package [T7] to generate a 
three-dimensional magnetic field map [5]. Granulation issues make the direct use of 
this map problematic in our particle stepper, so we use the map to find the parameters 
of an analytic model of the vector magnetic potential. A, from which we then derive the 
field. Using this analytic expression for A is computationally efficient, requires little 
memory, and eliminates the granulation issues. Over most of the particle-accessible 
space the fields derived from A are an excellent match to the numeric fields; the 
deviation between the numeric and analytic fields is never greater than about 2%, 
and is this large only near the axial ends of the octupole where particles rarely reach. 
However, while the fields derived from A satisfy V • B = exactly, they do not quite 
satisfy V x B = 0, and require the existence of unphysical currents, principally near 
the mirror coils. These currents are very small; over the majority of the trap, the 
unphysical current densities are more than four orders of magnitude lower than the 
typical current densities in the mirror coils. Even near the wall under the mirror 
coils where the unphysical current densities are largest, they are still more than two 
orders of magnitude lower than the typical current densities in the mirror coils. To 
further test the validity of this analytic calculation of B, we studied the distribution 
of annihilation locations with a computationally slower, but more accurate B found 
via the Biot-Savart line integral methodology. Since none of these studies showed 
statistically significant differences in the antiproton annihilation location distributions, 
we used the faster analytic calculation of B = V x A throughout this paper. Routines 
to calculate A and B were implemented independently in two different computer 
languages. The results of the two implementations were each checked against the 
numeric field map and against each other. The details of the calculation of A are 
given in [Appendix A[ 

An important advantage of the vector magnetic potential formulation is that it 
makes it trivial to calculate the electric field induced by the decaying magnetic field 
during the magnet shutdowns. This electric field, given by E = ~dA/dt, plays a 
key role in antiproton dynamics as it is responsible for conserving the third (area) 
adiabatic invariant |18] . 

The steady-state coil currents are measured to 1% accuracy, and this sets the 
accuracy to which the fields are known. During the magnet shutdown, the coil currents 
decay in a near exponential fashion with measured time constants near 9 ms. (In the 
simulations, we use the measured coil current decays to capture the small deviations 
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from exponential decay.) However, the changing magnetic field induces currents in 
the trap electrodes which retard the decay of the field. We have found these decay 
currents using the COMSOL Multiphysics package [TB] and a precise model of the 
electrode mechanical structure. The eddy currents depend on the resistivity of the 6082 
aluminum from which the electrodes are fabricated. This resistivity is 3.92 x 10~^ J7m 
at room temperature, and is reduced, at cryogenic temperatures, by the Residual 
Resistance Ratio (RRR), which we measured to be 3.06. We find that the eddy 
currents delay the decay of the magnetic field in a manner well-modeled by passing 
the coil-created magnetic field through a single-pole low pass filter; the filter time 
constants are 1.5 ms for the mirror coils and 0.15 ms for the octupole field. The eddy 
currents have more influence on the mirror fields than the octupole fields because the 
breaks between the electrodes do not interrupt the largely azimuthal currents induced 
by the mirrors, but do interrupt the largely axial currents induced by the octupole. 
The simulations employ these filters to model the effects of the eddy currents. 



2.2. Antiproton Simulations 

The antiproton simulations push particles in response to the Lorentz force ([2]) using 
the fields of Section |2.1[ Two codes were developed. The first and primary code 
propagates the full Lorentz force equations for the position r and velocity v using the 
Boris method 

r(i+|) = r(i-|)+<5M0 (7) 

v(i + 5i)=v(t)-^{E[r(t+|),<+|] 

v(t + St)+v(t) , St. St^. 

+ - ^ ^ xB[r(t+-),i+-]}, (8) 

where r is the antiproton position. This algorithm is an order St^ method for a single 
timestep St. It conserves the perpendicular energy exactly in a uniform, static field; 
this is particularly important as the simulations must conserve /ip adiabatically. 

The second code uses guiding center approximations, including E x B, curvature, 
and grad-B drifts, and propagates particles using an adaptive Runge-Kutta stepper. 
The results of the two codes were compared, and no significant differences were 
observed. Typical antiproton trajectories are described in [Appendix B| 

2.3. Antihydrogen Simulations 

The antihydrogen simulations pushed particles in response to ([T]) in the fields of 
Section |2.1[ Two adaptive Runge-Kutta stepper codes were developed independently 
and the results compared. No significant differences were observed. In addition, 
the usual convergence tests of the simulation results as a function of the time step 
were satisfactorily performed. Similar tests were also performed for the antiproton 
simulations. Typical antihydrogen atom trajectories are described in [Appendix C| 



3. Antiproton distribution and clearing 

Immediately after a mixing cycle, we axially "dump" the antiprotons and positrons 
onto beamstops where they annihilate. The dumps use a series of electric field pulses, 
and are designed to facilitate counting of the charged particles. They employ relatively 
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Figure 2. a) The total on-axis magnetic field B{z), and the electrostatic potential 
energy of an antiproton in the strong clearing fields at times tg and ti. b) The 
pseudopotential |9jl for antiprotons with perpendicular energy £±0 = 15, 24.4, 
and 50 eV. A well exists in the pseudopotential only for £±q > 24.4 eV. 



weak electric fields. (We switched from an "original" dump sequence to an "improved," 
more efhcient, dump sequence midway through the runs reported in this paper.) After 
the dumps, all the electrodes are grounded; any antiprotons that still remain in the 
trap must be trapped by the mirror and octupole fields alone. The magnitude of the 
mirror fields is plotted in Figure [2^. Next, we attempt to "clear" any such mirror- 
trapped antiprotons with a series of four clearing cycles. These clearing cycles use 
much larger electric fields than the dump pulses; there are two initial "weak" clears, 
and two final "strong" clears. The electrostatic potentials used in the strong clears 
are graphed as —qV{z,tQ) and —qV{z,ti) in Figure [2]i; the weak clear fields are half 
as large as the strong clear fields. 

A mirror-trapped antiproton can be thought to move in a pseudopotential $ 
which combines the electrostatic potential with an effective potential which derives 
from the invariance of , 

$(z,<) = -<zF(z,<)+/^pB(z). (9) 

For simplicity, we consider <i> on the r — axis only. Figure [2]d plots the 
pseudopotential for /ip/So = £±o = 15 eV, 24.4 eV, and 50 eV. For an antiproton 
to be trapped, a well must exist in the pseudopotential. This condition, which is a 
function of the perpendicular energy £±q only, replaces the prior trapping condition, 
^±o/^\\Q > ^pc in the presence of an electric field. For our parameters, a well only 
develops for antiprotons with £±o > 24.4 eV. Any antiproton with £±o < 24.4 eV 
will necessarily be expelled from the system by the strong clear field even if it has 
f||o = OeV. 

It might appear that antiprotons with £± > 24.4 eV would be trapped. But 
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Figure 3. The antiproton distributions that survive the clearing sequence 
simulations for the initial Distributions 1 and 2 defined in the text. Other initial 
distributions (not shown), which, for instance, start all the antiprotons close to 
the trap axis, yield similar thresholds. 



Figure [2|d shows the static pseudopotential; in the experiment, the clearing field swings 
from the potential shown in Figure [2^ at time to to the potential at time ti and back 
eight times (the first four swings, during the weak clears, are at half voltage). Each 
of these eight stages lasts 12 ms. Extensive computer simulation studies show that 
these swings expel all antiprotons with £± < £_LMirTrap = 50 eV. Two such studies are 
shown in Figure [3] 

The simulations are initiated with a postulated antiproton distribution before 
the clears. Unfortunately, we do not know this distribution experimentally (see 
Section[4|, so we use two trial distributions that cover the plausible possibilities: Both 
distributions assume a spatially uniform antiproton density throughout the trap region 
but differ in their velocity distribution. Distribution 1 has a velocity distribution that 
is isotropic and flat up to a total energy of 75 eV, while Distribution 2 has a velocity 
distribution that is isotropic and thermal with a temperature of 30 eV. Note that these 
distributions are intended to reveal the properties of antiprotons that could survive 
the clears. They are not intended to be representative of (and, in fact, are thought to 
be far more extreme than) the actual antiprotons in the trap. 

For both distributions, less than 2% of the antiprotons survive the clearing cycles 
and remain in the trap, and all that survive have £± > 50 eV. Further, those with 
£± > 50 eV are only trapped if they have very little fy, as not much parallel energy 
is needed for them to surmount the relatively shallow pseudopotential wells. For 
example, in Figure [2]d, the potential well for antiprotons with £± = 50 eV is only 
about 12 eV deep. 

The improved efficacy of the time-dependent clearing cycles over the static 
clearing potential comes from two factors: 1) The repeated voltage swings accelerate 
the antiprotons, in some cases non-adiabatically. This often gives them sufficient 
parallel energy to escape. 2) The potentials depicted in Figure [2^ are generated by 
voltages impressed on 21 electrodes. Four central electrodes have a significantly slower 
temporal response than the outer electrodes; this creates a momentary well that lifts 
and eventually dumps antiprotons with increased parallel energy, again raising the 
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Table 1. The average number of antiproton annihilations detected during the 
clearing operations. The data includes the false counts from cosmic background, 
which is separately measured, and given on the last row. The error is the statistical 
error of the average. The Trapping rows were measured during normal trapping 
attempts. The Benchmarking row was measured while deliberately creating high 
perpendicular energy antiprotons (see Sec.|5]l. The Full column shows the number 
of counts observed during the entire 24 ms time period taken by each clearing cycle. 
The Windowed column shows the number of counts between 0.6 ms and 2 ms in 
each cycle. We know from other data, not shown, that while trapping, almost all 
the antiprotons escape in this window. This is expected as it takes 2 ms for the 
clearing potentials to reach their peak. (Employing a 1.4 ms window increases the 
signal to noise ratio.) For the Benchmarking trials, antiprotons escape during the 
entire clearing cycle, and windowing would cut legitimate data. This data was 
collected by our detector in a non-imaging mode, wherein the detection efficiency 
is 70%-95% assuming most of the antiprotons hit near the trapping region. 



Full 



Windowed 



Trials 



Trapping 869 
-Original dumps 

First clear (Weak) 31.43±0.21 31.014±0.207 

Second clear (Weak) 0.38±0.02 0.022±0.005 

Third clear (Strong) 0.37±0.02 0.016±0.004 

Fourth clear (Strong) 0.31±0.02 0.022±0.005 

Trapping 371 
-Improved dumps 

First clear (Weak) 0.55±0.04 0.205±0.024 

Second clear (Weak) 0.34±0.03 0.035±0.010 

Third clear (Strong) 0.33±0.03 0.042±0.009 

Fourth clear (Strong) 0.24±0.03 O.OllitO.005 



Benchmarking 

First clear (Weak) 2460±150 

Second clear (Weak) 466±41 

Third clear (Strong) 283±30 

Fourth clear (Strong) 45.9±6.7 



27 



Background 



0.32±0.03 



0.019±0.002 



likelihood that they escape. 

We monitor the antiproton losses in our experiments during the clearing cycles 
(see table [l]) With the original dumps, a substantial number of antiprotons escape 
in the first clear. A few antiprotons escape during the second and third clears, but, 
to the statistical significance of the measurement, none escape in the last clear. With 
the improved dumps, far fewer escape in the first clear, a few, perhaps, in the second 
and third, and none in the last. It is telling that there is no upwards jump in the 
number that escape between the second and third clears (between the last weak and 
the first strong clear), as this lack suggests that there is no continuous distribution 
of antiprotons with a significant population with energies between £± w 25 eV, which 
are cleared by the weak clears alone, and £± = 50 eV, which are cleared by the strong 
clears. Thus, in conjunction with the simulations, we conclude that it is not likely 
that antiprotons with perpendicular energy less than 50 eV survive the clears and, 
therefore, none are likely present during the magnet shutdowns. 
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4. Mirror-trapped antiproton creation 

In this section we will describe three scenarios that could result in the creation 
of mirror-trapped antiprotons: creation during the initial capture and cooling of 
antiprotons from the AD; creation during the mixing of antiprotons into the positrons; 
and creation by the ionization of antihydrogen. We will show that none of these 
mechanisms are likely to produce mirror-trapped antiprotons with £± exceeding 50 eV. 
However, the calculations are sufficiently uncertain that they cannot guarantee that 
none are created. Instead, we rely on two other arguments: (1) As will be discussed 
in Section [7j the temporal-spatial characteristics of the candidate events are not 
compatible with mirror-trapped antiprotons. (2) By heating the positron plasma, 
we can shut off the production of antihydrogen |20j . When we do this, we observe 
essentially no trapped antihydrogen candidates (one candidate in 246 trials, as opposed 
to 38 candidates in 335 trials in [1]). The temperature to which we heat the positrons, 
approximately 0.1 eV, is negligible compared to the energy scales discussed in this 
section, and would have no effect on any mirror-trapped antiprotons created. These 
experiments are described in P" and will not be further discussed here. Taken together, 
these arguments allow us to conclude that few, if any, mirror-trapped antiprotons 
survive to the magnet shutdown stage where they could confound our antihydrogen 
signal. 

4.I. Creation on capture from the AD 

The AD delivers a short pulse of 5 MeV antiprotons to the ALPHA apparatus. 
These antiprotons are passed through a thin metal foil degrader, resulting in a broad 
antiproton energy distribution. The slowest of these antiprotons are then captured 
in a 3 T solenoidal field (eventually reduced to 1 T) by the fast manipulation of the 
potentials of a 3.4 kV electrostatic well [3T1 IHj. Once captured, about 50% of the 
antiprotons are cooled to several hundred Kelvin by collisions with the electrons in 
a pure-electron plasma that had been previously loaded into the same well [22]. The 
electrons themselves cool by emitting cyclotron radiation. The remaining 50% of the 
antiprotons do not cool: they are trapped on field lines at radii greater than the outer 
radius of the electron plasma and, thus, do not suffer collisions with the electrons. 
These uncooled antiprotons are removed from the trap by decreasing the trap depth 
to, ultimately, about 9 V on the trap axis, corresponding to 30 V at the trap wall. (The 
trap depth on the axis is less than at the wall because of the finite length to radius 
ratio, 20.05 mm/22.28 mm, of the trap electrodes.) As all of these preparatory steps 
occur before the neutral trapping fields are erected, any antiproton with f y exceeding 
9 (30) eV will escape before the neutral trap fields are erected and, thus, will not be 
mirror trapped. 

In principle, there is a remote possibility that a high perpendicular energy 
{£± > 50 eV) antiproton might be largely outside the electron plasma, so that it 
is not strongly cooled, but would have a parallel energy sufficiently low [< 9 (30) cV] 
that it could be trapped in the electrostatic well. Certainly, as shown by SRIM 
[13] calculations, a few antiprotons leave the degrader with such skewed energies. 
However, the antiprotons must surmount a ~ 50 V blocking barrier to enter the well. 
An antiproton could be mirror trapped only if (1) it possessed a high initial £± and 
a high enough £|| so that it could pass over this barrier, and then (2) underwent 
one or more collisions that reduced its £|| to less than 9 (30) eV while leaving its £± 
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above 50 eV. Such an evolution is unlikely to result from collisions with electrons 
as the antiproton-electron mass ratio requires multiple collisions to effect a significant 
change to the antiproton energy, and such a collision sequence does not favor a skewed 
distribution. Furthermore, we know experimentally that those antiprotons that do not 
cool quickly essentially never cool; lengthening the cooling time beyond some tens of 
seconds does not significantly increase the fraction of antiprotons that are cooled. 
Thus, it is unlikely that an antiproton would cool just enough to leave it in a state 
that could be mirror trapped, but not so much that it cools entirely. 

Alternatively, the collisions required to leave a mirror-trapped antiproton might 
be with another antiproton or with a background neutral gas molecule. The density 
of these necessarily high radius antiprotons is low, and, if they are to be mirror 
trapped, their perpendicular and hence total energy is high. The exact parameters 
to use in an antiproton-antiproton collision calculation are unknown, but, under any 
scenario, only a few antiproton-antiproton collisions will take place during the 80s 
cooling time. For example, for a plausible density of energetic antiprotons of about 
10'' cm~^, the probability that one 500 eV antiproton would suffer one collision in 80s 
is approximately 10~^. Furthermore, only a small fraction of these collisions would 
leave the antiprotons with the required skewed energy distribution. 

The neutral gas density can be estimated from the antiproton annihilation rate, 
and is on the order of lO^cm"'^, if, as is likely, the background gas in our cryogenic 
trap is H2 . While this yields an antiproton-neutral collision rate that is higher than 
that for antiproton-antiproton collisions, the collision rate calculated by extensions of 
the methods in |24j and |25j is on the order of a few tens of microHertz per antiproton, 
making it unlikely (few 10~^) that an individual antiproton will suffer a collision that 
will leave it with the energies required to be trapped. Individual antiprotons do not 
suffer multiple collisions with neutrals on the relevant time scale. 

4-2. Creation during mixing 

Antihydrogen is generated by mixing antiproton and positron plasmas after the neutral 
trapping fields are erected. By this point in the experimental cycle, the two species 
are cold; the antiprotons are at temperatures of less than two hundred Kelvin, and 
the positrons are at temperatures less than one hundred Kelvin p^. The expected 
number of antiprotons with an energy exceeding £_LMiiTiap in a thermalized plasma of 
particles and temperature T is A'^exp(— £_LMiiTrap/fcBr), where is Boltzmann's 
constant. This number is completely negligible for the relevant temperatures. The 
antiproton temperature would have to be approximately 200 times greater ( ^ 3 cV) 
for there to be an expected value of one or more antiprotons with energy greater than 
50 eV amongst the ^30,000 antiprotons present in one mixing cycle. Thus, there is 
no chance that thermalization of the initial antiproton plasma could produce mirror- 
trapped antiprotons. 

During the mixing cycle, the axial motion of the antiprotons is autoresonantly 
excited PSIITT] to ease them out of their electrostatic well and into the positron plasma 
(see Figure|4^-b); the antiprotons phase- lock to a weak, downward-frequency-sweeping 
oscillating potential applied to a nearby electrode. The autoresonant drive has a 
maximum potential drop of 0.05 V on the trap axis (0.1 V at the wall), and there are 
approximately A^a = 300 drive cycles. Naively, one might think that there are enough 
cycles that the drive could excite antiprotons up to the maximum confining potential 

* All gases but H, H2 and He freeze out; monatomic H is rare, and there is no source of He. 
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of 21 V on the trap wall. In reality, the antiprotons phase-lock at near 90° such 
that the impulse conveyed to the antiprotons on each cycle is small The typical 
antiproton gains just enough energy to enter the positrons: about 0.5 V on the trap axis 
when plasma self-fields are included. If, as occasionally happens, an antiproton loses 
phase-lock, it will gain a limited amount of additional energy stochastically in rough 
proportion to y/N]^. Further, this is axial energy; if the antiproton were to somehow 
gain more than 21 eV it would be lost immediately unless it had also experienced a 
sufficient number of collisions to posses substantial perpendicular energy. Under no 
scenario can the antiproton gain energy close to 50 eV of perpendicular energy directly 
from the autoresonant drive. 

The autoresonant process injects most of the antiprotons into the positrons, but 
some are left in the original side well with axial energies up to the electrostatic 
well depth of about 0.5 V near the trap axis. As mixing progresses, antiproton- 
antiproton collisions cause additional antiprotons to fall into this side well and into 
the electrostatic well on the other side of the positron plasma (see Figure |4j;) . As 
there is no direct mechanism to transport these antiprotons radially outward |29) . 
most will remain at or near their original radius (between 0.4 and 0.8 mm depending 
on the details of the procedures in use at the time.) Approximately 50% of the 
particles eventually fall into the two side wells, so the number of antiprotons in the 
side wells eventually approaches the un-mixed antiproton number. Measurements on 
similar plasmas show that they thermalize in times on the order of the one second 
that the mixing continues [22 . (Unlike in Section 4.1 the density of these near-axis 
antiprotons is relatively high.) Measurements also show that evaporative cooling will 
set their temperature to several times less than the well depth [TT]. Thus, the near-axis 
antiproton temperature in the side wells is considerably less than 0.5 eV. The expected 
number of antiprotons in such a plasma having a perpendicular energy greater than 
50 eV is negligible. 



4-3. Creation by ionization of antihydrogen 

Antihydrogen in the ALPHA experiment is believed formed largely by three body 
recombination. This process creates the atoms in highly excited states that can be 
ionized by sufficiently strong electric fields [TJl [IS] . The strongest electric fields in our 
trap are found close to the trap wall at the electrode boundaries, and can be as large 
as ii^max = 42 V/mmj5] A newly ionized antiproton will be accelerated by these fields, 
and can pick up perpendicular energy. However, a careful map of the electric and 
magnetic fields over the entire trap shows that the perpendicular energy gain cannot 
exceed more than 3eV before the antiproton settles into its E x B motion, so this 
process cannot lead to mirror-trapped antiprotons. 

The arguments in the two previous subsections strongly suggest that antihydrogen 
cannot be born with substantial center-of-mass kinetic energy under our experimental 
conditions. If an antihydrogen atom were, nonetheless, somehow born with high 
kinetic energy, this energy would be conveyed to the antiproton upon ionization. 
Naively, this could lead to a mirror-trapped antiproton. However, there is an upper 
limit to the amount of energy an antiproton could possess after ionization. The 
limit comes from the Lorentz force equation ([2]). A particle moving at velocity v 

tt Very close to the electrode gaps, the electrode corners will increase the field beyond -Emax = 
42V/mm. However, any antiproton born close enough to corners to feel this enhancement will 
almost surely hit the wall immediately. 
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Figure 4. On-axis electrostatic potentials in the mixing region of our apparatus. 
The green dots are a cartoon depiction of the evolution of the antiprotons. a) 
Before the autoresonant injection of the antiprotons. Note how the positron space 
charge flattens the vacuum potential, b) Immediately after autoresonant injection 
of the antiprotons. c) At the end of the mixing cycle. 



perpendicular to a magnetic field B feels a force that is equivalent to that from an 
electric field of magnitude vB. This magnetic force qvB can ionize an antiproton just 
as can an electric force qE. Thus, if an antihydrogen atom is sufficiently excited that 
it can be ionized by the large electric field of strength i?max or less near the trap 
wall, it will always be ionized by passage through the magnetic field at the center 
of the trap where it is created if it is moving faster than approximately Enmx/Bo- 
This sets a rough upper limit on the maximum kinetic energy that a high radius, 
newly ionized antiproton can have of less than 10 eV. If an antihydrogen atom has 
more kinetic energy, it will either: 1) be in a relatively low excited state such that it 
will not be ionized at all, and will hit the trap wall and annihilate promptly. Or (2) 
be in an ionizable state and be ionized close to the trap axis by the magnetic force, 
where it will be thermalized and cooled by the abundant population of antiprotons 
and positrons found there. A more exact calculation, given in [Appendix D lowers 
this bound substantially for most antiatoms. 

The side wells near the trap wall are as deep as 21V. An antiproton that fell into 
one of these side wells, either indirectly by ionization or directly by some unknown 
process during mixing, could pick up substantial parallel energy. However, the 
density of antiprotons is very low at large radii, and antiproton-antiproton collisions 
proportionally infrequent. Multiple collisions would be required to transform the 
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Figure 5. The distribution of tlio deliberately-created mirror-trapped 
antiprotons, shown just after the antiprotons were dropped over the 40 V potential 
barrier, at 120 s post drop, and when thermalization was ended at 420 s post 
drop. These distributions were obtained by slowly lowering one of the confining 
electrostatic barriers, and measuring the number of escaping antiprotons as a 
function of the barrier height. 

maximum parallel energy of 21 eV into perpendicular energy of more than 50 eV. 
Collisions with neutrals, of course, can only lower the antiproton energy. Thus, we 
can conclude the parallel energy possessed by an antiproton cannot be converted into 
sufficient perpendicular energy to lead to mirror trapping. 

5. Antiproton simulation benchmarking during magnet shutdowns 

The arguments in the previous sections suggest that there are few, if any, mirror- 
trapped antiprotons. This tentative conclusion relies on information gleaned from the 
simulations of the efficacy of the clearing cycles. Ultimately, though, we rule out the 
existence of mirror-trapped antiprotons by comparing their simulated post magnet 
shutdown dynamics to our experimental observations. A direct, independent test of 
the simulations powerfully buttresses our conclusions. 

We performed such a test by deliberately creating a population of mirror-trapped 
antiprotons. We began by capturing approximately 70,000 antiprotons from the AD. 
These antiprotons were injected over a potential barrier into a deep well, giving them 
a parallel energy £j| of approximately 40 eV. Then the antiprotons were held for 
420 s, during which time collisions partially thermalized the populations, transferring 
parallel energy into perpendicular energy £±. The mean antiproton orbit radius also 
expanded during this time to approximately 1.5 mm. Initial, intermediate, and final 
f II distributions are shown in Figure [5j We have no independent measure of the £± 
distribution. 

After the thermalization period, the octupole and mirror coils were energized, 
followed by the removal of the electrostatic well that had been confining the 
antiprotons. Once this well was removed, the antiprotons remaining in the system 
must have been mirror trapped. However, many of these antiprotons were not deeply 
mirror trapped < £'_LMirTrap ~ 50 eV), and, as can be seen in "Benchmarking" 
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Figure 6. Comparison of tlie z-t anniliilation locations of mirror-trapped 
antiprotons (symbols) and the antiproton simulations (dots), a) Compares "Left 
Bias" (blue upward pointing triangles and blue dots), "No Bias" (green circles 
and green dots) , and "Right Bias" (red downward pointing triangles and red dots) 
for normal current decay times, b) Compares the "No Bias" dataset in a) with 
normal shutdown timing, to a No Bias dataset in which the octupole decay onset 
was slowed by about 7 ms (purple squares and purple dots) relative to the mirror 
decay onset. The annihilations near z = —183 mm, and ±137 mm are at radial 
steps in the trap wall. The detector resolution was approximately 5 mm in z, and 
100 lis in t; the simulation points were randomly smeared by these resolutions. 



grouping m table [1] many were expelled during the clears. 

After the clears, the magnets were turned off, and the annihilation times t and 
positions z of the remaining antiprotons were recorded. The results of 27 of these cycles 
are shown in Figure[6^. During most of these cycles, a Bias electric field (see Figure]?]) 
was applied during the magnet shutdown whose intent was to aid the discrimination 
between bare antiprotons and antihydrogen; the charged antiprotons should be pushed 
by the Bias field so that they preferentially annihilate on the right side (Right Bias) or 
on the left side (Left Bias) of the trap, while the uncharged antihydrogen atoms should 
be unaffected by the Bias field. In addition, the Bias fields make the pseudopotential 
wells shallower, so the antiprotons escape and annihilate sooner than when no bias is 
applied. Figure [HJd shows the effect of delaying the octupole shutdown onset by about 
7ms relative to the mirror shutdown onset. 

Also plotted in Figure [6] are the results of simulating 3364 post-clear survivors. 
Since we can only characterize the pre-clear and magnet shutdown antiprotons 
imperfectly (see Figure [s]), we must make an estimate for the distribution to use 
in the simulation. We believe Distribution 2, defined in Section [3j is most appropriate 
as it has a plausible temperature and no strict upper bound on £± . Figure [6] shows 
that simulations match the experimental data well. Thus, we can confidently use 
the antiproton simulations as a tool to aid in the discrimination between mirror- 



Antihydrogen and mirror-trapped antiproton discrimination 



17 



300 




-200 -100 100 200 

z (mm) 



Figure 7. The electrostatic Bias potentials —qV(z), and the on-axis 
pseudopotentials ^(z) for £±^ = 60 eV for full and half strength mirror fields 
and for the a) Left Bias, b) No Bias, c) Right Bias cases. When the Bias fields 
are applied, the antiprotons are localized at the ends of the trap. The localization 
is preserved as the magnets lose strength during the magnet shutdown. 



trapped antiprotons and antihydrogen. These tests also confirm that the bias fields 
work as expected. The antihydrogen simulation uses the same magnetic field model 
as the antiproton simulation, so we have benchmarked the field component of the 
antihydrogen simulation as well. 

6. Postulated Antihydrogen Energy Distribution 

As described in Section |4.2[ antihydrogen atoms are created by mixing antiprotons 
with positrons. Initially, the antiprotons have more kinetic energy than the positrons, 
but we estimate that the antiprotons come into thermal equilibrium with the positrons 
before the recombination occurs. The positron density is 5 x 10'' cm^'^ and the positron 
temperature is 40K [2 . We use [SIT to compute a slowing rate of ^ 200 s~^. From |31| . 
the three body recombination rate is approximately 0.1 s^^, but this is the steady state 
rate to reach a binding of 8k^T. Because antihydrogen atoms that have a binding 
energy of Ik^T will mostly survive the fields of our trap, the recombination rate 
will be approximately ten times higher. This is in approximate agreement with our 
measurements. Consequently, we expect that the antiprotons cool to the positron 
temperature before forming antihydrogen. 

Because the positron mass is negligible compared to the antiproton mass, a 
newly formed antihydrogen atom inherits its centcr-of-mass kinetic energy from the 
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antiproton from which it is formed. Thus, we expect that the antihydrogen itself 
is in thermal equilibrium with the positrons, and possesses the same distribution 
function — except that the trapped antihydrogen distribution function is truncated at 
the energy of the neutral trap depth, fxrap = 0.54 K. The positron temperature 
is much greater than this energy. Consequently, we expect that the velocity space 
distribution function f{v) is essentially flat over the relevant energy range for the 
trapped antihydrogen atoms, and the number of atoms in some velocity range dv is 
f{v)v'^dv (X dw oc d£ . The number of atoms trapped should be proportional 
to ^^xrap- Note that because f{v) is essentially flat in the relevant region, the 
antihydrogen distribution v'^dv, once normalized, does not depend in any significant 
way on the temperature of the positrons. However, for concreteness, we did our 
principal antihydrogen simulations with a temperature of 54 K. 

The simulations reveal that the energy distribution is not strictly truncated at 
the trapping depth (see Figure [8^) [5]. There exist "quasi-trapped" stable trajectories 
with energies up to about 0.65 K; similar trajectories exist in neutron traps |32) . 
Quasi-trapped trajectories exist because the antiatom motion is three dimensional. 
Rarely is all of the antiatom's motion parallel to the gradient of |B| at the orbit 
reflection points at high |B|. Any motion perpendicular to V|B|, and the kinetic 
energy associated with this perpendicular motion, is not available to help penetrate 
through the reflection point. Hence, the antiatom may be confined even if its energy 
exceeds the maximum trapping depth. Being only quasi-trapped, these antiatoms are 
more susceptible to perturbations than antiatoms trapped below the trapping depth. 
We do not know if the quasi-trapped trajectories are long-term stable. 

The positron plasma is Maxwellian in the frame that rotates with the positron 
plasma. This rotation modifies the lab frame distribution. If the positron density 
were to be very high, the rotation would impart significant additional kinetic energy 
to the antiprotons, and, hence to the resulting antiatoms. This would result in fewer 
antiatoms being caught in the trap. For our densities and fields, however, this effect is 
small. The reduction in the number of antiatoms that can be trapped from this effect 
is less than 5%. 

7. Trapping Experiments 

During the 2010 experimental campaign, we observed 309 annihilation events 
compatible with trapped antihydrogen. These events were observed under a number 
of different conditions, including runs with Left Bias, No Bias, and Right Bias, and 
with the antihydrogen held for times ranging from 172 ms to 2000 s. The conditions 
under which the observed events were obtained are listed in tabled 

Figure |9] plots the spatial and temporal (z-t) locations of the observed 
annihilations after the octupole and mirror fields were turned off. The figure compares 
the observed annihilation locations to the locations predicted by the antihydrogen and 
antiproton simulations. The initial distributions in these simulations were the flat 
antihydrogen distribution discussed in Section |6j and the antiproton Distribution 1 
defined in Section |3] We chose Distribution 1 here rather than Distribution 2, because 
we wanted to maximize the number of antiprotons just above the mirror-trapping 
barrier £_LMirTrap; for the benchmarking test in Section [5j we chose Distribution 2 
because we had independent evidence (Figure [5]) of the existence of antiprotons well 
above fj^MiiTrap- However, as is evident from comparing Figure [6^ and Figure [9ji, the 
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Figure 8. a) The antihydrogen energy distribution f{£). The solid green area is 
a histogram of the energy of the trapped antihydrogen atoms as predicted by the 
simulation from a starting population of atoms at 54 K. The blue line plots the 
expected E^^"^ dependence up to an energy of 0.54 K. This line ends at the vertical 
gray line, past which point all the antihydrogen atoms are quasi-trapped. The 
red points plot the energy distribution function reconstructed from the observed 
data. The reconstruction process is discussed in |Appcndix E| The error bars come 
from Monte Carlo simulations of the reconstruction process and represent only 
the statistical errors, b) The time of annihilation after the magnet shutdown as a 
function of the initial energy for simulated antihydrogen atoms. (For clarity, only 
a representative 2000 point sample of the 35000 simulated antiatoms is plotted.) 
The function of the gray band is described in [Appendix E[ c) Histograms of the 
number of annihilations as a function of time after the magnet shutdown, as 
observed in the experiment (red points) and in the simulation (solid green area.) 
The error bars on the experimental points come from counting statistics. 



differences between the annihilation locations for these two distributions are minor; 
principally, some of the higher energy antiprotons in Distribution 2 annihilate closer 
to the center of the trap than the antiprotons in Distribution 1. 

In general, the agreement between the observed events and the antihydrogen 
simulations is excellent; in contrast, the vast majority of observed events are 
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Table 2. Observed trapping events during the 2010 experimental campaign. 
"Hold Time" denotes the time interval between when most of the antiprotons 
were dumped from the trap, thereby ceasing antihydrogen synthesis, and when 
the trap magnets are turned off; i.e. the approximate minimum time that the 
antihydrogen was trapped. As the trapping rate improved continuously during 
the experimental campaign in 2010, and long Hold trials were all clustered near 
the end of the campaign, no conclusions about the lifetime of antihydrogen in 
our trap can be reached from from the ratio of observed trapping events to the 
number of trials [2]. 



Hold 
Time (s) 


Left Bias 


No Bias 


Right Bias 


Total 


Trials 


0.2 


73 


41 


13 


127 


613 


0.4 


129 




17 


146 


264 


10.4 


6 






6 


6 


50.4 


4 






4 


13 


180.4 


10 




4 


14 


32 


600.4 






4 


4 


38 


1000.4 


5 




2 


7 


16 


2000.4 






1 


1 


3 


3600.4 










1 


Total 
Trials 


227 
577 


41 
227 


41 
182 


309 


985 



incompatible with the antiproton simulations. As expected, the locations of the 
observed events are independent of the bias electric fields, as they are in the 
antihydrogen simulations. The simulations show, however, that the annihilation 
locations of postulated antiprotons are strongly dependent on the bias fields. Other 
simulations, not shown here, show that these conclusions remain true in the face of 
antiproton energies up to several keV and gross magnetic field errors. 

In [Tl[51|3] we limited our analysis to annihilations which occur within 30 ms of the 
beginning of the magnet shutdown. This criterion was based on the observation that, 
in the simulations, 99% of the antihydrogen atoms annihilated by 30 ms. Here, and 
in P], we impose an additional requirement that \z\ < 200mm; beyond this region, 
the efficiency of the detector falls and the accuracy of the detector reconstructions 
becomes suspect (see Figure |9]:,f). In the first 50 ms after the shutdown, we observed 
four events which fail these cuts, and, thus, do not appear to be antihydrogen atoms. 
These events also appear to be incompatible with the antiproton simulations. While 
we have no definitive explanation of these events, there are several possibilities: 1) 
Even if all 309 events were genuinely due to antihydrogen, we would expect 1% of the 
events to be improperly excluded because of the t criterion; the total number of events 
thereby improperly excluded would be expected to be three. 2) As discussed in [Tl[3], 
cosmic rays are miscategorized as antiproton annihilations at a rate of approximately 
47 mHz. The events discussed here were observed in approximately 985 x 50 ms « 50 s, 
so we would expect to observe approximately 2 such miscategorized cosmic rays, some 
of which could occur outside the cut boundary. 3) The basic z resolution of our 
detector is approximately 5 mm, but there is a low probability long tail of badly 
resolved annihilations (Figure [9j^) . Some of these observed events may be outside 
the \z\ < 200mm window because they were poorly resolved. 4) The trap electrodes 
have offsets of up to about 50 mV due to the non- ideal behaviour of the electrode 
amplifiers. This creates shallow wells, which might store antiprotons outside of the 
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Figure 9. a) Spatial and temporal {z-t) locations of the annihilations during the 
trapping events, and the annihilation coordinates predicted by the antihydrogen 
simulations (small gray dots). Table [2] details the trapping conditions, b) z 
histograms of the annihilation locations. The observed locations agree well with 
the predictions of the antihydrogen simulation, and are independent of the Bias 
conditions, c) Detector efficiency as a function of z, as calculated by GEANT 3 
1331 . d) Similar to a), but with the annihilations predicted by the antiproton 
simulations for Left Bias conditions (left clump of purple dots), No Bias conditions 
(central clump of green dots) and Right Bias conditions (right clump of red dots), 
e) Similar to b), but with histograms from the antiproton simulations. The counts 
in the simulation histograms are divided by a factor of five so that the observed 
event histogram is also visible.) f) Percentage of the reconstructions that are more 
than 50 mm and 100 mm from their true position, as calculated by GEANT 3. 

region in which the clearing fields are applied, and which might cause antiprotons to 
be released at odd times and positions. 

As remarked above, annihilations typically occur within 30 ms of the magnet 
shutdown. The time history of these annihilations contains information about the 
energy distribution of the antihydrogen atoms |2] • Figure [sjj plots the annihilation 
time as a function of energy as found in the simulations. As expected, the higher energy 
antiatoms, which are freed at higher values of the diminishing trap depth, annihilate 
sooner than low energy antiatoms. Figure ^ shows a histogram of the expected and 
observed annihilation times. The observed points are well predicted by the simulations. 
From the data in Figure [8];, the original energy distribution of the antiatoms can 
be coarsely reconstructed, as shown in Figure [8^. To within the predictive power 
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of the reconstruction, the energy distribution follows the expected £^^^ plus quasi- 
trapped distribution. The reconstruction algorithm, and its very significant limits, are 
described in |Appendix E[ The influence of the energy distribution on the z distribution 
is described in 



8. Conclusions 



We have presented a detailed study of the behaviour of antihydrogen atoms and 
antiprotons confined in a magnetic minimum trap. This study was used to 
guide experiments which eliminate antiprotons as a possible background in recent 
antihydrogen trapping experiments. We have demonstrated how the very different 
behaviour of the neutral and charged particles lead to very different loss patterns in 
time and space when the magnetic minimum trap is rapidly de-energized. These 
different loss patterns have been a crucial factor in the identification of trapped 
antihydrogen. Finally, we have shown how we can use the simulations to reconstruct 
the energy distribution of the trapped antihydrogen from the time history of the loss 
after de-energization. These studies and tools have provided important insights into 
the nature of antihydrogen trapping dynamics. 

In the future, it may be possible to discriminate between antihydrogen and 
antiprotons via a resonant interaction with the atomic structure of the antiatoms. 
Such resonant interactions could photoionize the antiatoms or flip their spins such 
that the antiatoms become high field seekers. The technique of field-ionization, which 
has been successfully used with excited antiatoms [H], does not work with ground 
state antiatoms because the fields required to strip a ground state antihydrogen atom 
are too large, and, thus, would not detect the long-trapped atoms discussed here [5]. 
Until efficient resonant interactions with the antiatoms can be obtained, the techniques 
demonstrated in this paper will remain a crucial tool in the endeavour to increase the 
trapping rates and pursue the path towards detailed spectroscopy of antihydrogen. 

This work was supported by CNPq, FINEP/REN AFAE (Brazfl), ISF (Israel), 
MEXT (Japan), FNU (Denmark), VR (Sweden), NSERC, NRC/TRIUMF, AITF, 
FQRNT (Canada), DOE, NSF (USA), and EPSRC, the Royal Society and the 
Leverhulme Trust (UK). 



Appendix A. Magnetic Field Formulas 

In this Appendix, we develop the analytic model of the magnetic field referred to in 
Section |2.1| and required for use with the simulations. 



Appendix A.l. Mirror coils 

By comparison with the precise Biot-Savart fields, we found that the magnetic field 
from each individual mirror coil could be accurately approximated using a pair of 
circular loops. From Jackson [33], the vector potential from a single loop is 

^^^^ (r2 + a2)3/2 (1 + y (^2 + a2)2 + (^-1) 

where 4> = arctan(?;/x), C — IfiQC? /A: is a constant, a is the radius of the loop, 
— a;;2 _|_ g^j^j-j ^2 _ ,.2 _^ ^2 .^yj^j^ x,y,z measured from the center of the circle 
defined by the loop. Unfortunately, the series converges very slowly near the mirror 
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and this formula, although accurate, was abandoned. Instead, we used a method based 
on guessing a form for A. The guess is inspired by the form of the exact A from a 
single loop: 

A^ ^ C-^[{a^ + r^- 2aAs)-i/2 _ (^2 ^ ^2 ^ 2a\s)-^/% (A.2) 

where all of the parameters are as before and A is a dimensionless fit parameter. 
Note, the choice A — \/3/2 — 0.866 exactly reproduces the first two terms of the 



exact A^ (A.ll for a single loop. The two mirrors are slightly different. Our fit gave 
a = 45.238 mm, A = 0.9019 and a loop separation of 8.251mm between the two coils 
of the left mirror, and A = 0.9027 and a loop separation of 8.579 mm between the two 
coils for the right mirror, and a separation between the two mirrors of 274 mm. We 
found that these choices gave max(|Bfit — BcxactI) < 0.02 T when the mirror field was 
1 T. This maximum error occurred on the wall of the trap directly underneath the 
mirrors; for \/x^~+^ < 15mm the maximum error was ^ 0.01 T. 

Appendix A.2. Octupole field 

The vector potential for an infinite octupole is 

Aoo = Fs^ cos(40)z, (A.3) 

where s = \/ ^ ip- and _F is a constant. A finite, symmetric octupole can be written 
as 

A, = {Fi{z)s^ + F6(z)s« + ^^8(z)/ + ...) cos(40), (A.4) 
where the i^'s are functions to be determined later. The condition 

V'^A^ = (A.5) 
gives the relations 

Fe= - — , (A.6) 
fa 20' ^ ' 

^« - - 4I - W' ^^-'^ 

etc. 

In order to satisfy V • A = there must be non-zero components of A in the s 
and <p directions: 

As - (G5(z)s^ + G7(z)s^ + ...) cos(40), (A.8) 
A^ = (i/5(2)s^ + Hr{z)s^ + ...) sin(40). (A.9) 
The G's and i/'s are determined by the equations 

Id,,, 1 dAj, dA^ , , , 

V ■K=Q=- — {sAs) + -^ + —^, A.IO 
s OS s ocp oz 

V^A^ = 0, (A.ll) 

V'^Ay = 0. (A.12) 

The second two relations lead to = H^, Gt = Hj etc. The first relation leads to 

G,{z)= -^Fiiz), (A.13) 

Gr(z) = - —F'(z) = (A.14) 
7V ^ 12 ^ 240 

etc. 
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Note that there is only one free function, ^4(2); aU of the other functions are derivatives 
of this one. To get a fit to F4 we need a function that looks like a symmetric plateau. 
We chose to use the complementary error function 

Fi{z) = D[erfc((z - Zf)/Az) - erfc((z + Zf)/Az)], (A.15) 

where D is a, constant, ±Zf are the approximate ends of the octupole and Az is the 
distance in z over which the octupole drops to ^ 0. 

In our fit, we found Zf — 129.46 mm and Az = 16.449 mm. This form was able 
to get max(|Bfit — BexactI) < 0.02 T when the field at the wall was ^ 1.5 T. Because 
of the way in which the functions were chosen, the condition V • A is always exactly 
satisfied if the orders of the expansion are kept the same in all three components of A. 
The condition V^A = is satisfied only to the extent that enough terms are retained 
in the expansion. For our parameters, V^A is small in the region of interest inside 
the trap. 



Appendix B. Mirror- Trapped Antiproton Trajectories 

Electrostatically trapped antiprotons follow regular trajectories similar to those shown 
in Figure |B1[ The antiprotons oscillate between the two ends of the electrostatic well, 
following field lines that typically extend between a radial minimum at one end of the 
well to a radial maximum at the other end of the well. These radial maxima occur 
in magnetic cusps [351 four to each side, caused by the octupole's radial fields. 
Guiding center drifts cause the antiprotons to slowly rotate around the trap axis, so 
that the trajectories slowly alternate between cusps at each end. The consequences 
of this motion, like the existence of a limit on the maximum allowed well length, have 
been explored in a series of papers [32l |38l ESI |39l [36] . 

Mirror-trapped antiprotons trace far more complicated trajectories, as shown in 
Figure |B2| Typically, the z motion follows a relatively slow macro-oscillation that 
extends over the full axial extent, and a faster micro-oscillation, over a more limited 
axial extent. Each micro-oscillation typically travels between two large-radius, local 
octupole cusps, though sometimes the micro-oscillation ends at a low-radius point 
near one of the mirror coils. Since the only mechanism for reversing the antiproton 
motion is a magnetic mirror reflection, the reversal necessarily occurs at a relatively 
large value of |B|. Indeed, the reflection always occurs at the same value of |B|: at 
the field magnitude at which all of the antiproton's kinetic energy is completely tied 
up in its conserved magnetic moment fip (see Figure |B3]and Q) 



As can also be seen in Figure |B3[ the trajectories take the antiprotons closest 
to the trap wall in the center of the trap. The electric field sloshing in the clearing 
cycles leaves the antiprotons with on the order of 5-10 eV. Consequently, the 
antiprotons oscillate from one end of the trap to the other rapidly; for the trajectory 



in Figure B2 the macro-oscillation bounce frequency is on the order of 20 kHz. After 
the magnet shutdown, antiprotons escape over a time of more than 10 ms; thus, the 
antiprotons typically make hundreds of bounces during the shutdown process. This 
allows the antiprotons to find the "hole" in the trap center, and causes the antiproton 
annihilations to be concentrated there in the No Bias case (see Figure |6]) . When a 
Bias is applied, the center of the pseudopotential moves to the side (Figure [?]), and 
the annihilation center follows. 
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Figure Bl. Typical trajectory of an electrostatically trapped antiproton, with 
£|| = £± = 10 eV and a starting radius of 11mm. The magnetic fields are held 
constant at the values given in Section ^ The antiproton oscillates in a well 
formed by two end electrodes biased to —140 V, separated by a 80.2 mm grounded 
electrode. As in Figure^ the trap's central axis points along z, and the center of 
the trap, at 2; = 0, is in the center of the grounded electrode, a) Axial {t—z), b) 
radial (z-r), and c) transverse (x—y) projections of the motion. 



Appendix C. Minimum-B Trapped Antihydrogen Trajectories 



A typical minimum-B trapped antihydrogen atom trajectory is graphed in Figure CI 
In the transverse plane, the antiatom oscillates radially, with a varying rotational 
velocity; a Fourier transform (not shown) of the x or y motion yields a broad range of 
frequencies. This is expected as an r^ potential in the transverse plane, to which the 
potential in our trap approaches, is known to yield chaotic motion [40;. 

Typically, the antiatom trajectories cover the transverse plane reasonably 
uniformly, with little azimuthal structure, but are peaked at the outer radial edge 



where the antihydrogen atoms reflect (see Figure C2). The ultimate goal of these 
experiments is, of course, to use spectroscopy to search for differences between 
antihydrogen and normal hydrogen. The plots in Figure |C2| suggest that the 
trajectories do not sharply constrain the waist of a probe laser or microwave beam. 

The axial motion is quasi-harmonic with a well defined oscillation frequency that 
typically remains constant for many oscillations. Occasionally, as shown in Figure [C3j 
the frequency jumps due to interactions with the transverse motion. Since the z- 
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Figure B2. Typical trajectory of a mirror-trapped antiproton, with £"11 = 10 eV, 
£± = 60 eV and a starting radius of 11mm. All electrodes are grounded, a), b) 
and c) are described in Figure [bT] 



oscillation frequency is on the order of 100 Hz, the antiatoms bounce only a few times 
during the magnet shutdown. Unlike mirror-trapped antiprotons, the antiatoms do 
not have time to find the low |B| hole in the 2;-center of the trap, and, consequently, 
they annihilate over a broad region in z. 

Appendix D. Ionization of fast antihydrogen 

To study the ionization of fast moving antihydrogen atoms, such antiatoms were 
propagated in a constant axial magnetic field B — IT, and a radial electric field 
E = neqp/2€o = EqP arising from the space charge of the positron plasma. Here, 
He = 5.5 X 10^ cm~^ is the plasma density, and p — {x, y, 0). The equations of motion 
in terms of the center of mass coordinates Rcm: V^Mi a-nd the relative coordinates r, 
V, are 

M'VcM = qv^B + qEap (D.l) 

= g(VcM + Av) X B + 9£;o(pcM + Ap) +Fe (D.2) 

= (?Av X B + qEeff + Fc, 

where M is the total mass of the atom, /i the reduced mass, A = {mp — me)/M and 
Fc the Coulomb force. The effective electric field Eeff = Vcm x B -|- qEq {pcM + Ap) 
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Figure B3. The same trajectory as in Figure [B2]3, but plotted for 1ms rather 
than 0.1ms. The solid and dashed lines are lines of constant |B|, plotted at the 
two angles, 22.5° and —22.5°, of the octupole cusps. One set of lines is plotted at 
|B| = 1.26 T, the reflection field for the plotted trajectory. The other set is plotted 
at the value of |B| such that the largest radial extent equals the wall radius Ry/ 
at 10 ms after the magnet shutdown, a typical time for an antiproton to hit the 
wall. The inset figure shows that the trajectories terminate on one angle or the 
other depending on their z direction. 



is the sum of the regular electric field and a term proportional to the center-of-mass 
velocity of the atom. 

The coupled equations (D.ll and (D.2) were solved using an adaptive step size 
Runge-Kutta algorithm. The antihydrogen atoms were initialized at a trap radius of 
0.5 mm and some given initial binding energy S^, and initial kinetic energy fkin in the 
transverse plane. Binding energy is here defined as in the field-free situation, i.e. as 
the sum of the kinetic energy of the positron and the Coulomb potential. For each 
parameter set {fbji^^kin}, one thousand trajectories were calculated. Each trajectory 
was followed for a maximum of 2 fis or until the atom was ionized. The fraction of 
trajectories leading to ionization, as well as the time until ionization, were recorded. 

The magnetic field creates an effective harmonic confinement for the positron in 
the transverse plane. Hence, strictly speaking, one cannot have field-ionization (in the 
sense that r — >■ oo), unless there is also some axial electric field present, which was 
not the case in these simulations. However, a strong radial electric field will induce 
a positron-antiproton separation much larger than the atomic size in the field-free 
situation. Such a positron will be bound only by a negligible binding energy, and the 
antiatom will almost instantly be destroyed by either collisions with another positron 
(inside the positron plasma) or by a weak axial electric field (just outside the plasma). 
We regard any antihydrogen atom bound by less than 2 K, which is much less than 
the plasma temperature, as ionized. 

The fraction of antihydrogen trajectories leading to ionization is shown in 
Figure |D1| for various initial binding energies and center-of-mass velocities. An 
antihydrogen atom is stable against ionization by an axial electric field for binding 
energies £b > 2\/ (g^/47reo)g£^z- Typical electric fields in the trap are of order 10 V/cm, 
corresponding to stability for f b ^ 30 K. Any antihydrogen atom with a binding 
energy less than 30 K will be ionized by the effective electric field with more than 
99% efficiency at even moderate kinetic energies of 0.1 eV. However, very close to the 
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Figure CI. Typical trajectory of a minimum-B trapped antihydrogen atom. 
The antiatom started with a kinetic energy of 0.5 K. a), b) and c) are described 
in Figure [BT] 
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Figure C2. Transverse projections of three typical antihydrogen atom 
trajectories, each for antiatoms with 0.5 K energy, but with differing, randomly 
picked, initial directions. The antiatoms were propagated for 100 s. Each 
projection was scaled to the same maximum on the linear colourmap. Below 
each projection is the corresponding density profile. 

electrode boundaries, the electric fields can be much larger, corresponding to stability 
only for > 150 K. Our simulations show that such deeply bound atoms will require 
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Figure Dl. Probability of ionization witiiin 2 /js as a function of kinetic energy 
for antiiiydrogcn atoms witii the listed binding energies £^,. 



much larger kinetic energies to ionize in the lower-field region in the center of the trap 
(see Table Im]). 
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Table Dl. Minimum kinetic energy ^Tj^in of an antihydrogen atom required 
for ionization within 2^s with 90% probabihty (column 2) and 99% probability 
(column 3) for different binding energies ^Tb. 
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Appendix E. Energy Reconstruction 

The trapped antihydrogen energy distribution function /(£) can be crudely 
reconstructed from the time history of the annihilations after the magnet shutdown. As 
shown in Figure [Sja, antiatoms of a given energy £ annihilate over a broad distribution 
of times. The overall probability distribution function for the antiatoms to annihilate 
at time t can be found by integrating the probability P{t\£) of annihilation at time t 
of antiatoms with specific energy £ over the antiatom energy distribution function: 

fit)= r A£P{t\£)f{£). (E.l) 

This equation can be exploited by guessing a distribution function /(f), calculating 
P{t\£) with simulations, and comparing it to a histogram of the observed data. This 
"forward" method was explored in [5] and in Figure [s]:. 
Alternately, we can write 

/•oo 

f{£)= / dtP{£\t)f{t). (E.2) 
In this appendix we explore the consequences of employing this "inverse" equation. 



We perform the integral in (E.2 1 as follows: (a) For each annihilation event, construct 
a narrow band around the annihilation time in the antihydrogen simulation results. A 
typical such band is shown in gray in Figure [8]d. (b) From this band, randomly select 
a fixed number of the simulated annihilation events. (We selected 20 such samples 
in the reconstructions in this paper.) This effectively finds and samples P{£\t). (c) 
Aggregate all the energies from the randomly selected samples for each observed event, 
effectively integrating over t as properly weighted by f{t). (d) Construct the histogram 
of these aggregated samples; this is the reconstructed energy distribution. 

The eight subgraphs of Figure |E1| show a study of the reconstruction process 
for eight trial distributions. For each trial distribution, we analyzed two thousand 
sets of 309 Monte-Carlo generated annihilation events, each event obeying the trial 
distribution particular to the figure subgraph. The average over all of the resulting 
reconstructions for each subgraph is then plotted. One can see that the reconstruction 



is coarse. Figures El i-e show that the mean energy of the distribution is recovered 



approximately, as well as some features of the higher moments of the distribution, but. 
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for more pathological distributions, Figures |El| -h show that the ability to recover 
these higher moments is limited. 

The reconstruction errors stem from two causes: (1) The band of energies at 
each time is broad (see Figure This results in sharp features being smeared; 



this problem is particularly relevant in Figures El '-h. (2) The reconstruction has 



a difRcult-to-quantify memory of the original distribution used in the simulations 
underlying the reconstruction; this problem is particularly acute in Figures [El| 3-c. The 
reconstruction can be improved by iteration; an initial reconstruction, done employing 
the original simulation results, can be used to determine the approximate temperature 
of the experimental data, and then the reconstruction rerun using a simulation with 
a more appropriate temperature. 
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Figure El. Montc-Carlo study of the trapped antihydrogen distribution function 
reconstruction algorithm. In all eight subgraphs, the green, solid region is a 
histogram of a postulated distribution function f{€) whose reconstruction (red, 
blue and pink lines) is being attempted. The histograms are generated from: 
a) the distribution of surviving (i.e. trapped) antihydrogen atoms as predicted 
by the antihydrogen simulation from a starting population of atoms at 54 K. 
This is the distribution principally studied in this paper. In this and in all 
subsequent cases, the histograms are not smooth because they are generated from 
a finite number of samples from the starting population; b) the distribution of 
surviving antihydrogen atoms starting from a population of atoms at O.IK; c) 
the distribution of surviving antihydrogen atoms starting from a population of 
atoms at 0.01 K; d) a distribution similar to that in a), but with the distribution 
artificially forced to be flat out to the trapping energy, and then rolled off with the 
same quasi-bound distribution as found in a); e) a distribution that is similar to 
d), but which is artificially forced to increase linearly out to the trapping energy; 
f) the quasi-bound antiatoms in a) only; g) the non-quasi-bound antiatoms in a) 
only, h) a double humped distribution. In a), the red dotted line is the average 
reconstructed distribution function found using an inversion based on the 54 K 
simulation study, as described in [Appendix E[ For comparative purposes, this red 
dotted line is replicated in all the subgraphs. In all the subgraphs but the first, the 
red solid line is the reconstruction of the postulated distribution function in the 
particular subgraph, also found with an inversion based on the 54 K simulation 
study. In b) and c), the blue dashed lines are reconstructions found with an 
inversion based on a 100 mK simulation study. In c), the pink dashed line is the 
reconstruction found with an inversion based on a 10 mK simulation study. In 
each of the plots, the reconstructions are averaged over two thousand 309 point 
Monte-Carlo generated event sets. (Figure shows the typical reconstruction 
and error with just one 309 point event set — our actual data.) From this survey, 
it is clear that the mean energy of the distribution, as well as some coarse features 
of the distribution, are recovered, but sharp features are lost. 
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